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NUMERICAL SIMULATION OF THE SOFIA FLOWFIELD 


Abstract 

This report provides a concise summary of the contribution of computational 
fluid dynamics (CFD) to the SOFIA (Stratospheric Observatory for Infrared 
Astronomy) project at NASA Ames and presents results obtained from closed- 
and open-cavity SOFIA simulations. The aircraft platform is a Boeing 747SP 
and these are the first SOFIA simulations run with the aircraft empennage 
included in the geometry database. In the open-cavity runs the telescope is 
mounted behind the wings. Results suggest that the cavity markedly influences 
the mean pressure distribution on empennage surfaces and that 110-140 
decibel (db) sound pressure levels are typical in the cavity and on the horizontal 
and vertical stabilizers. A strong source of sound was found to exist on the rim of 
the open telescope cavity. The presence of this source suggests that additional 
design work needs to be performed in order to minimize the sound emanating 
from that location. A fluid dynamic analysis of the engine plumes is also 
contained in this report. The analysis was part of an effort to quantify the 
degradation of telescope performance resulting from the proximity of the port 
engine exhaust plumes to the open telescope bay. 
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plume spread angle 
ratio of specific heats, 1 .4 
plume spread rate 
normalized radial dimension 
velocity ratio 
atmospheric density 
plume spread parameter 


fluctuating (turbulent) value 


hot engine exhaust 

potential core of the fan plume 

chord 

centerline 

potential core of the hot engine exhaust plume 
fan 

fan to hot engine exhaust (ratio) 
full-width-at-half-maximum 
free-stream to fan exhaust (ratio) 
maximum 

single-equivalent jet 
thermal 

aircraft coordinate system (x is measured along the 
longitudinal axis of the fuselage from the nose of the 
aircraft to the tail, y is measured from the fuselage 
out the starboard wing, and z points in the vertical 
direction) 

telescope coordinate system; and z'" 

pass through the center of the telescope air bearing 
and are parallel to the aircraft coordinates x, y, and z, 
respectively 
free-stream 

location of the half-maximum 


Introduction 

This report describes computational fluid dynamics (CFD) research in the 
SOFIA (Stratospheric Observatory For Infrared Astronomy) project at the NASA 
Ames Research Center. The research has entailed using CFD technology for 
simulating the flow field in the vicinity of SOFIA. Computations ultimately will be 
used to assess effects of the open telescope cavity on the fuel consumption, in- 
flight safety, and fatigue of the aircraft platform (a Boeing 747). An important part 
of the work completed to date has involved validating CFD capability in 
transonic cavity-flow applications by comparing computational results to wind- 
tunnel experiments. CFD has also been used to evaluate whether mounting the 
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"P e in a cavity behind the wings is technically viable, rather than forward 
of the wings and behind the cockpit. The latter was the first configuration 
proposed for SOFIA. a 


This report contains results of open and closed-cavity SOFIA computations. The 
closed-cavity and one of the open-cavity simulations were run at 747SP cruise 
flight conditions. A third simulation was run at conditions of wind-tunnel 
experiments performed during 1994 at NASA Ames Research Center. The 
telescope is mounted behind the wings and the empennage is included in the 
geometry database. As noted in subsequent discussion, the impact of the open 

♦ cavity . on ^ 0W ’ n vicin 'ty of the horizontal and vertical 

stabilizers is of importance to design engineers concerned with control of the 
aircraft platform and fatigue of empennage surfaces. 


Background 

Since it began five years ago the SOFIA CFD program has been tasked with 
accomplishing four objectives. The first entailed developing and validating CFD 
capability in transonic cavity-flow applications and has been largely 
The second involved simulating the flow field in the vicinity of 
oUrlA in order to assess effects of the open telescope cavity on the sensitivity 
and resolution of the infrared telescope (i.e., its so-called seeing) and on aircraft 
performance (i.e., its control, safety, and fuel consumption). The first part of this 
objective has been largely met as well. The third objective involves addressinq 
technical issues regarding effects of the open cavity on fatigue of the aircraft 
a * n u.-it” e( i tS u of di ^ erent cavity-door designs on telescope performance and 
stability of the cavity-port shear layer. This objective has not been accomplished 
and will remain the focus of future effort. The last objective follows immediately 
from the first three. Once it becomes operational, SOFIA will be based at NASA 
Ames. During the projected twenty-year lifetime of the SOFIA program a 
validated design tool will be needed to assess effects of design changes which 
might be contemplated for improving telescope seeing or for enhancing aircraft 
performance. Since wind-tunnel experiments are much more costly and time- 

consuming to run than computer simulations, CFD is the logical choice for that 
design tool. 


SOFIA was one of fifteen projects chosen in 1991 by the National Research 
Council as crucial to continued progress in astronomical science during the 
next decade. SOFIA is planned as a successor to the Kuiper Airborne 
Observatory (KAO) currently based at Ames and will consist of an infrared, 2.5 
!T#u ter Cassegrain telescope mounted in a Boeing 747SP or Boeing 747-200 
When it becomes operational, the capability of SOFIA will be significantly 
greater than the KAO's and SOFIA will offer advantages which earth- and 
space-based instruments cannot. Specifically, SOFIA will operate at an altitude 
“Jove m ?,f t . of the ab sorbing water vapor in the atmosphere. Consequently, 
SOFIA will be capable of collecting data which is not measurable from the 
earth s surface. Furthermore, SOFIA will function as a mobile observatory which 
can be used to collect data from transient astronomical events at remote 
locations around the world. Although a space-based telescope would 
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undoubtedly possess more sensitivity, SOFIA’S telescope will cost considerably 
less and have a larger diameter mirror with greater resolution. Moreover, SOFIA 
will be considerably less expensive to service, since maintaining, upgrading, 
and modifying its instrumentation will be comparatively easier, less costly, and 
less time-consuming. As has the design of the KAO, the design of SOFIA will 
allow installation of state-of-the-art instruments throughout the operating life of 
the aircraft, and will provide a test platform for new instrumentation which may 
eventually be used on space-based platforms. 

Nevertheless, the use of an aircraft as an observation platform creates a 
number of technically challenging problems. Solid materials typically have a 
very narrow transmittance bandwidth for infrared radiation. Consequently, since 
the spectrum of radiation from the sky and celestial sources is very broad, the 
SOFIA telescope will be mounted in a fuselage cavity which will be uncovered 
and open to the atmosphere when data is being collected. The effect of the 
open cavity on the in-flight performance and safety of the aircraft, as well as on 
the optical resolution and pointing accuracy of the instrument, were technical 
questions which arose in the preliminary design of SOFIA. Fortunately, the 
design criteria for optimizing aircraft performance are the same as those for 
optimizing telescope performance. Performance of the telescope will be 
enhanced by a cavity configuration which results in a stable, attached cavity- 
port shear layer. Aircraft fatigue will also be reduced by such a design, and 
control of the aircraft will be improved as well. 

The issue of an aft telescope mount arose in 1991 as part of an effort to identify 
a less expensive, yet technically viable aircraft cavity and telescope design. In 
the SOFIA design first proposed, the telescope was mounted forward of the 
wings immediately behind the aircraft cockpit. However, mounting the telescope 
behind the wings has distinct cost advantages because an aft mount will require 
a simpler cavity door design, installation of one less pressure bulkhead in the 
aircraft, and fewer modifications of aircraft control and power systems. 
Nevertheless, technical problems which are not severe in a forward-mount 
configuration may be important and difficult to overcome with a telescope 
mounted farther aft. The dynamic steadiness of the cavity-port shear layer and 
acoustic levels in the cavity are influenced by the thickness of the fuselage 
boundary layer upstream of the cavity aperture. Since the boundary layer is 
much thicker and more turbulent behind the wings than forward of them, the 
shear layer across the cavity port is also thicker and more turbulent in an aft- 
mount configuration. Consequently, shear-layer turbulence levels (as measured 
by turbulent kinetic energy) will be greater and light passing through the thicker 
layer will be diffracted more. Most of the SOFIA design effort completed to date 
has presumed that a Boeing 747SP would be the aircraft platform for the SOFIA 
telescope. However, if the telescope is mounted behind the wings, the SP may 
not be suitable for this purpose. There are two reasons. First, an adverse 
pressure gradient exists on the SP fuselage immediately downstream of the 
cavity and a region of separated flow is evident in the aft-cavity CFD simulations 
run to date. Second, the cavity aperture on an SP cannot be fitted with an 
external door because movement of the door is restricted by the aircraft's 
vertical stabilizer. An SP can be fitted with an internal cavity door, but an 
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internal door restricts space in the cavity and decreases the allowable focal 
length of the telescope. These issues have recently led the SOFIA project to 
consider the Boeing 747-200 as a possible platform for SOFIA. 

To date CFD has been an integral part of a general test and analysis program in 
the SOFIA project. CFD has provided a unique tool for evaluating design 
concepts and has provided estimates of unsteady telescope loads, which were 
not measured in the 1990 wind-tunnel experiments conducted at NASA Ames. 
A second set of experiments was completed during the summer of 1994. The 
focus of these experiments was the aft-mount telescope configuration and a 
comparison of the flow field characteristics in the vicinity of the telescope cavity 
mounted on both the 747SP and 747-200 aircraft platforms. 

Wind-tunnel experiments were run in 1990 on a sting-mounted, clipped-wing 
aircraft model and were used to validate NASA Ames CFD capabilities in the 
SOFIA project. The experiments focused on investigating resonance 
characteristics of the open telescope bay in a forward-mount configuration and 
were tasked specifically with obtaining a design which results in a stable shear 
layer across the cavity port at transonic Mach numbers. A stable shear layer 
reduces density fluctuations in the optical path of the telescope and aero- 
acoustic and aerodynamic fatigue resulting from pressure fluctuations in and 
around the open telescope bay. Moreover, effects of flow-field oscillations on 
the empennage and, consequently, control of the aircraft are undoubtedly 
mitigated when the shear layer is stable. A cavity-port ramp design which 
generates a stable shear layer was identified in the experiments and numerical 
simulations captured the essential characteristics of the flow-field resonance 
generated by that design, thereby validating the CFD effort. The ramp will 
probably need to be modified when the telescope is mounted behind the wings, 
in part because the direction of flow over the cavity will be different in that 
configuration. Moreover, the installation of an aperture door will undoubtedly 
impact the design of the ramp. 

The experiments in 1990 measured only surface pressures in the telescope 
cavity, but CFD can be used to examine the physics of the flow in a more 
detailed manner than experiments can. Indeed, numerical simulations provide 
information which experiments generally cannot - integrated pressure data, 
from which aerodynamic loads on the telescope can be computed, temperature 
and density data in the cavity and across the cavity-port shear layer, by which 
thermal and density effects on the aero-optics of the SOFIA telescope may be 
examined, and fluctuating pressures in and around the telescope cavity, by 
which aero-acoustic fatigue of the aircraft may be quantified. (CFD has already 
provided estimates of unsteady telescope loads, which were not measured in 
the experiments.) CFD can therefore be expected to have a major impact on the 
future direction and progress of the SOFIA program, since experiments are 
costly and cannot be planned or scheduled without significant lead times. 

The O-type cavity-port design identified in the 1990 experiments has been used 
in all subsequent numerical and experimental work in the SOFIA project. 
However, an important result of the open-cavity simulation described in this 
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report points to the need for design modifications in order to eliminate a strong 
source of acoustic energy on the rim of the open telescope bay. 

Geometry 

Boeing provided NASA with a complete 747SP surface definition in the fall of 
1992. The data included the geometry of the aircraft empennage, which had not 
been available before. All of the CFD simulations completed prior to that time 
were run with a sting-mounted fuselage and geometry derived from sources 
secondary to Boeing, i.e., a CAD definition of the Microcraft wind-tunnel model 
used in the 1990 experiments, a coarse 747-200 surface definition which has 
been used in panel code simulations, United Technologies Pratt and Whitney 
JT9D engine sketches, and NASA blueprints. The Boeing data was used to 
verify location, length, and angle-of-attack of the wings, the definition of the 
wing-root fairings, and fuselage taper in the region behind the wings and in 
front of the empennage where the cavity is to be located in the aft-mount 
configuration. Data was not used to refine the surface definition of the wingtips 
and nacelle and engine geometries, which were flow-through in the data bases 
previously supplied to NASA, but which were modified for powered conditions 
from information in Boeing and Pratt and Whitney sketches. These parts of the 
surface definition are not of critical interest to the design of the cavity. 
Unfortunately, the surface definitions provided by Boeing included only 
unloaded wing geometry. As a result, the loaded wing geometry of the panel 
code definition was used in the simulations described in this report, but no 
attempt was made to verify that the loading was appropriate for 747SP cruise 
conditions. 

AtwoodT2,3,4 use d GRIDGEN2D and HYPGEN software to construct the surface 
and volume grids in his SOFIA simulations and most of these grids did not need 
to be modified for use in the CFD simulations. Atwood subsequently constructed 
empennage grids from the Boeing data obtained in 1992 and used them to start 
a closed-cavity computation. This simulation was not completed, but its solution 
was used to initialize the flow field in the closed-cavity flight simulation 
described in this report. 

All of the CFD simulations were run on a 7% model of the aircraft platform. The 
empennage surfaces on the starboard side of the aircraft were generated by 
reflecting Atwood's half-model grids about the y = 0 plane. The telescope 
elevation angle, measured from the horizontal, was 40°. Before the wind-tunnel 
simulation was run all engines and nacelles were removed from the geometry 
database and the wings were clipped outboard of the position where the 
outboard engines would normally be mounted. This configuration was the same 
as that tested in the wind-tunnel experiments run during 1994. 

Method 

Two simulations described in this report were run at SOFIA cruise conditions (M 
= 0 85 at an altitude of 41,000 feet) and one at conditions of the wind-tunnel 
experiments. In the former the engines were powered and simulations included 
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the engine inflow and outflow boundary conditions used by Atwood 3 , for which 
j^u an< ^ Wh,tne y P rovided engine thrusts and mass flow rates through the fan 
and hot engine cores. The Pratt and Whitney data were also used to compute 
tail cone surface temperatures. 


All simulations were run on the Cray C-90 supercomputers in the Central 
£?/cDr. m ™? nd Numer,cal Aerodynamic Simulation Facilities at NASA Ames 
OVERFLOW was used to compute solutions. OVERFLOW contains implicit' 
approximately-factored algorithms and has been used to run all SOFIA cavitv- 
flow simulations completed to date. The diagonalized Beam-Warminq 
numerical scheme and the Baldwin-Lomax turbulence model in OVERFLOW 
were used in all calculations. The effects of turbulence in the cavity-port shear 
layer were simulated by the Baldwin-Lomax shear-layer model, which has been 
described by Atwood 3 . 


The closed-cavity simulation was run on a computational domain containing 
thirty-eight grids and approximately 1.1 x 10 3 grid points. The domain in the 
open-cavity simulation consisted of forty-five grids and approximately 2 5 x 10 6 
grid points. An overset grid interface routine (PEGASUS) was used to establish 
mte Jl"9 ri ° communication and compute the grid-boundary interpolation 
coefficients required by OVERFLOW. H 

The open-cavity simulations were run time-accurately. The closed-cavity 
solution provided initial conditions at all points in the computational domain of 
ttie open-cavity flight simulation, with the exception of the flow field in the cavitv 
The cavity flow field was initialized by Atwood's 3 - 4 simulations, which were run 
on cavity grids identical to those in this computation. The time step was 
approximately 1.35 x 10' 4 seconds in this run. A data record one(1) second 
long on which oscillations with frequencies less than two(2) hertz cannot be 
resolved, requires approximately 7.4 x 10 3 time steps. Since the CPU cost to 
run OVERFLOW is approximately 7 x 10" 6 sec per time step per grid point a 
one-second data record requires approximately thirty-five C-90 CPU hours to 
compute solutions on domains containing 2.5 x 10* grid points. The data record 
of the open-cavity wind-tunnel simulation was approximately 0.074 seconds 
long. A 0.074 record length at wind-tunnel conditions is equivalent to a 1 2 
second record at cruise flight conditions. 

If they exist, oscillations with frequencies of 0(1) hertz (at flight conditions) are of 
particular concern to design engineers in the SOFIA project. Telescope motions 
excited by flow field phenomena at these frequencies can be eliminated bv 
proper design of the tracking mechanisms. However, a one-second data record 
is too short for computing reliable statistics of motions at these frequencies. 

Results 

Closed-Cavity Simulation, Flight Conditions 

The simulations reported by Atwood 3 - 4 were used to initialize the closed-cavity 
run, which was the first SOFIA CFD computation to include the aircraft 
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empennage in the geometry database. The simulation was run approximately 
3.8 x 10 3 steps after initialization until the C p distributions on the fuselage had 
achieved steady-state equilibrium. Fig. 1 contains plots of experimental and 
computed, steady-state C p distributions on the top, port side, and bottom o e 
fuselage during flight. The experimental data was provided by Boeing 6 and 
used by Atwood 3 - 4 to validate CFD results on a sting-mounted aircraft without 
empennage. However, the data was probably taken f r °m me^urements an 
altitude lower than the 41 ,000 foot cruise altitude of the SOFIA CFD simulations. 

For the purpose of plotting, the starboard side C p distribution is identical to that 
on the port side of the aircraft and is not shown in Fig. 1. However C p 
distributions on empennage stabilizer surfaces were slightly asymmetrical and 
suqqest that the simulation had not reached steady-state everywhere in the 
computational domain. Nevertheless, it is clear that CFD reproduces very well 
the pressure distribution on the fuselage and, in particular, in the empennage 
region, where one would expect the flow field to be strongly influenced by an 
open telescope cavity. 

Plume Fluid Dynamics 

If the telescope is mounted behind the wings, infrared radiation (IR) from the hot 
exhaust plumes of the port engines may scatter onto the primary mirror of the 
instrument and degrade its performance. This problem was the subject of an 
independent investigation by the SOFIA project, which was tasked with 
quantifying the degradation, if any. Part of the mve^igatior i involved 
experiments with the Shuttle Carrier Aircraft, which is a modified Boeing 747- 
100 powered by Pratt and Whitney JT9D engines - the engines most likely to be 
used to power SOFIA as well. Results of the investigation suggested that the IR 
emission from engine exhausts would not significantly degrade telescope 
performance, except possibly at very low telescope elevation angles (Dinger et 

al. 6 ). 

Attention in this report will focus on the exhaust plume of the port outboard 
engine, because thermal IR scattering from this engine to the telescope mirror is 
siqnificantly greater than that from the inboard engine. In the experiments the 
exhaust plumes of both inboard and outboard engines were Imaged by IR 
cameras mounted in the NASA Ames Lear Jet, which flew alongside the Shuttle 

Carrier (Dinger et al. 6 ). 

The brightness distribution of the jet exhaust is dictated by the mixing and 
subsequent cooling of the gases in the plume behind the aircraft Empirica^ 
correlations suggested by Antonia and Bilged, Morris 6 , Lau 9 , and Tanna and 
Morris 10 were used to construct a simple model describing the evolution of the 
plumes, ie„ their spread rate, length of potential core, and mean and turbulent 
temperature distributions. Model predictions were then compared with 
measurements of peak temperatures, plume full-width-at-half-maximum df wh m . 
and an estimate of the length of the hot potential core. M ° re /e cent 'y- r ® s c ulta 
from the closed-cavity CFD simulation were also compared with predictions o 
the model. All of the comparisons serve to validate CFD and the empirical 
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model. The validation provides a rational basis for their future use should 
questions again arise with regards to effects of SOFIA'S engine emissions on 
performance of the telescope. 


Appendix A contains a summary of the equations used in the empirical 
formulation of the model and Figs. 2 and 3 contain sketches of the model's 
construct. Mixing of the fan and hot exhausts along the tail cone was neglected 
and. the jet at the downstream tip of the cone was assumed to consist of a hot, 
circular core of area equal to the area of the nozzle exit, surrounded by a 
concentrically circular jet of area equivalent to the exhaust area of the fan. The 
radii r c and r f of the concentric disks are depicted in Fig. 2. Velocities and 
temperatures at this location (x' = 0) were assumed to be equal to the 
corresponding velocities Vf and V c , and temperatures 7> and T c , of the fan and 
engine exhausts, respectively, for which data was available from the 
experiments or was provided by Pratt and Whitney. This construct provided an 
upstream boundary condition for the model. 

Downstream of the potential core the plume evolves as a fully-developed jet 
with radial temperature and velocity distributions which are approximately 
Gaussian, but which can be represented by a sinusoidal distribution suggested 
by Morris . The thrust of the fully-developed jet was required to be the same as 
the thrust of the coannular jet at the x' = 0 upstream boundary (Tanna and 
Morris 10 ). This requirement leads to the concept of the so-called single- 
equivalent jet (Ko and Kwan 11 ), which provided an estimate of the cross- 
sectional area of the jet at the end of the potential core. The spread rate 
correlations in Appendix A cannot be used to estimate the size of the fully- 
developed jet at the end of the potential core, because the region between x' f 
and x cpt where the fan and hot engine exhaust plumes mix, is not self- 
preserving. The spread rate correlations strictly apply to self-preserving jets, 


Potential core lengths 

The radial spreading of the plume and the erosion of its constant-velocity and 
constant-temperature inner region are results of the intense mixing of the fan 
and hot engine exhaust gases, which occurs over considerable distance 
downstream of the engine nozzle. The centerline temperature decreases 
downstream of the potential core, whose length is a function of Mach number in 
the core gases and fan and engine temperatures, velocities, and sizes The 
model makes use of a correlation suggested by Tanna and Morris 10 to estimate 
the lengths x' f and x' p of the fan and hot exhaust potential cores. 

A potential core length x^ of 7.7 meters is predicted by the model for the 
nominal flight conditions of the Shuttle Carrier experiments, most of which were 
conducted at M = 0.6 and an altitude of 35,000 feet. (Some data were 
collected at an altitude of 26,000 feet, but no comparisons were made with 
these measurements.) This value corresponds well with the 7-8 meter length 
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apparent in the peak short-wave band (SWB) and long-wave band (LWB) 
temperature distributions reported by Dinger et al. 6 . 

Plume width 

In the region of the potential core, which ends at x / r c = 17.1, approximately 
(Fia 4) the model predicts a nearly-constant plume width dfwhm of 
approximately 0.9 meters for the conditions of the Shuttle Carrier experiments. 
This result agrees reasonably well with the LWB average of 0.7 meters, which is 
also virtually constant in the same region. Farther downstream the plume 
spreads at an approximate angle of 7°, but spreading is not apparent in the 
experimental data. 

CFD predictions of plume width from the closed-cavity simulation (at M = 0.85 
and an altitude of 41 ,000 feet) were also compared to the measurements. 
Results are plotted in Fig. 5. The difference in flight conditions has little effect on 
the results and agreement between the measurements and CFD predictions is 

good. 

Fig. 6 contains model and CFD predictions of plume width at SOFIA cruise 
conditions. According to the model, the potential core ends at x'/r c = 18.0 
( x ' _ 8.1 meters). Moreover, the model overpredicts the width of the plume, 
particularly in the fully-developed region downstream of this location. 
Nevertheless, the growth rate of the plume (i.e., the rate of change of the 
plume's width with respect to streamwise distance x') predicted by CFD is 
similar to that predicted by the model in this regime. 


Peak inferred temperature 

The peak temperature in the plume occurs along the jet centerline at any 
downstream location. In the potential core the temperature is constant and 
equal to 322°C, the static temperature of the hot engine exhaust. LWB 
measurements are in good agreement with this prediction, but SWB 
measurements are about 100°C less. The centerline temperature decreases as 
distance increases downstream of the potential core, but the decrease in the 
data is much more rapid. For example, a centerline (peak) temperature of 
300°C is predicted at a downstream distance of 12 meters, whereas the LWB 
data suggests that the peak temperature is less than 100°C at the same 
location. 


Peak rms temperature fluctuations 

Estimates of rms temperature fluctuations in the plume can be deduced from the 
well-known Crocco relationship (Cebeci and Smith’*) and the experimental 
data of Antonia and Bilger 7 , Lau 9 , and Smith and Hughes 13 . Appendix A 
contains an expression for the maximum rms temperature fluctuations in a tuliy- 
developed heated jet. The temperature fluctuations are expressed as a function 
of the velocity fluctuations, which are maximum at z' ~ z 0 .s- According to 
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experiment 7 - 13 , {Jv*) ma x - 0.3{V d -V„) at this location. Moreover, at SOFIA 

cruise conditions, f T - 1. Consequently, (Vr^) ma * = 0.3(T c/ -TJ. No attempt 
was made to deduce temperature fluctuations from the image data, although it 

is possible to do so. Nevertheless, an estimate of (Vt' 2 ) max , from which density 
fluctuations can be determined, is important in analysis of telescope seeing. 

Plume centerline velocity and temperature distributions 

Figs. 7 and 8 contain model and CFD predictions of velocity and temperature on 
the centerline of the exhaust plume at SOFIA cruise conditions. The centerline 
values represent maxima in the plume at any streamwise location x'. The model 
assumes that centerline values of velocity and temperature in the potential core 
are essentially undiminished from the values in the x' = 0 plane. However, at 
locations farther downstream, the model predicts a much faster rate of 
centerline decay 

Radial velocity and temperature distributions 

Figs. 9 through 14 contain plots of radial velocity and temperature distributions 
in the fully-developed region of the plume at SOFIA cruise conditions. 
Downwash from the wings distorts the CFD distributions, which would otherwise 
be symmetric about the z'= 0 axis. The CFD predictions shown were computed 
in the vertical x'z' plane. 

Tail Cone Temperatures 

All of the SWB and most LWB images in the experimental data were saturated 
in the vicinity of the engine tail cone downstream of the exhaust nozzle. As a 
result, the infrared emission of the cone and plume in this region could not be 
determined directly from most of the data. However, the brightness of the tail 
cone is a function of its surface temperature, which can be estimated using 
empirical heat transfer correlations. It is necessary only to assume 
thermodynamic equilibrium before computing the convective and radiative heat 
transfer between the cone and hot exhaust gas flowing over it. 

Radiative heat transfer rates were calculated assuming an emissivity of 1.0 for 
the cone and 0.1 for the exhaust gas. Empirical flat-plate correlations suggested 
by White 14 were used to compute skin friction and corresponding convective 
heat transfer coefficients. Reynolds number in the calculation of skin friction was 
based on distance from the turbine exit plane upstream of the engine exhaust 
nozzle. It is possible to correct coefficients computed from flat-plate correlations 
for the axisymmetric geometry of the tail cone. The correction involves 
computing an effective Reynolds number based on experiment (White 14 ). 
However, the correction has negligible effect on surface temperatures in these 
calculations and was not implemented. It is also possible to account for effects 
of compressibility (Mach number) on convection, but this correction is 
insignificant at the conditions of the Shuttle Carrier experiments. 



Table I contains the results of calculations for the outboard engine on the port 
side of the aircraft. Event numbers correspond to changes in engine throttle 
settings in the experiments. The surface temperatures are significantly higher 
than those obtained from LWB images and plotted in Fig. V-6 of Dinger et al. 6 . 
Moreover, surface temperatures deduced from the images are significantly 
lower (by 80°C, approximately) than the static temperatures of the hot exhaust 
gas at the exit of the engine nozzle. This result seems implausible, in part 
because the core temperature of the exhaust remains undiminished near the jet 
centerline well downstream of the tail cone. 

There is no reason to believe that the discrepancy between measured and 
computed surface temperatures is the result of poor assumptions or 
inaccuracies in the heat transfer calculations. For example, neglecting radiation 
emitted by the exhaust gas has little effect on computed temperatures, since the 
exhaust gas radiation is a small part of the energy budget. Moreover, a gas 
emissivity greater than 0.1 results in higher computed temperatures and a 
larger discrepancy between experimental and calculated values. 

If radiation or convection were entirely neglected in the calculations, surface 
temperatures would be well approximated by the adiabatic wall temperature in 
the different events. The latter temperatures are listed in Table I and are greater 
than the surface temperatures computed. However, radiative heat transfer is not 
negligible at the conditions of the Shuttle Carrier experiments and, in the 
absence of convection, computed and measured temperatures match only if the 
exhaust gas emissivity is greater than 1.0. This requirement is physically 
unrealistic. 

Nevertheless, empirical relationships for convective heat transfer coefficients 
often correlate poorly with available data and it is necessary to consider the 
possibility that inaccuracies in the calculation of this coefficient are responsible 
for the difference between measured and calculated tail cone temperatures. 
The temperatures deduced from the LWB images are correct only if convective 
heat transfer coefficients are approximately one order of magnitude smaller 
than the computed ones. Coefficients of such small magnitude are not 
consistent, however, with available heat transfer data. Consequently, one must 
conclude that the surface temperatures measured in the Shuttle Carrier 
experiments are inaccurate. 

Open-Cavity Simulation, Flight Conditions 

This simulation was initialized from the closed-cavity run on all grids in the 
computional domain, with the exception of grids in the cavity and in the 
immediate vicinity of the cavity aperture. The open-cavity run of Atwood (1993), 
which was computed on cavity grids identical to those in this simulation, 
provided the initial flow field for the cavity and cavity aperture grids. After 
initialization the simulation was run time-accurately for approximately 1.5 
seconds (real time) before data was collected. At SOFIA cruise conditions, a 
747 SP travels approximately seven(7) of its lengths in this period of time. This 


W ^ S deemed sufficient for generating a state of statistical quasi- 
equilibrium in the unsteady flow field in and around the open cavity and on 
empennage surfaces. The simulation was run another 1.1 seconds of real time 
dunng which data was collected for analysis. The 1.1 -second data record is 
sufficiently long for studying the essential characteristics of the unsteady flow 
Put is too short for calculating accurate statistics of flow field phenomena in a 
frequency range of a few hertz. Unsteady phenomena in this range are of 
particular importance to the SOFIA project. They will influence the design of the 
telescope air bearing and tracking mechanisms and impact estimates of fatigue 
loading when the cavity is open to the atmosphere. 

Figs. 15 through 22 contain sketches of the aircraft and show locations of the 
telescope cavity, cross-sections through the cavity and empennage surfaces 
a . nd P™ts wh ® re spectra of sound pressure levels have been computed. The 
athwartships plane AA through the fuselage and cavity in Figs. 15 and 16 is 
s own in Fig^ 20. Plane HH through the fore-and-aft symmetry plane of the 
tetescope in Fig. 20 is shown in Fig. 21 . Fig. 22 depicts the location of points on 
the forward and aft cavity bulkheads. The cross-sectional outline of the 
telescope shown in this figure serves to orient the view (cf. Fig. 20). 

Empennage C p Distributions 

The C p distributions in Figs. 23 through 25 were computed from 1.1 -second 

?hS! a - 9 ^ S ° f P r f s ® ures at Points on empennage surfaces. The distributions were 
then interpolated onto cross-sections parallel to the longitudinal axis of the 
fuselage. The interpolation was necessary because cross-sections in the 
computational domain are not generally oriented parallel to the fuselage axis 
Cross-sechons on the vertical stabilizer are nearly parallel to the axis, but C n 
distributions were interpolated on that surface as well, even thouqfi 
interpolation was not as important a consideration. 

^ C f c C % di ?L" bution . s in planes BB and cc throu 9 h the vertical stabilizer 
(Figs. 1? and 18) are shown in Fig. 23. Not surprisingly, there is asymmetry in 
the distributions when the telescope cavity is open. In plane BB the influence of 
the open cavity is predominately at mid-chord and farther aft. The influence in 
p ane CC is farther forward. The asymmetry in the distributions suggests that the 
open cavity affects the flow field over the entire span, and that effects are largely 
undiminished at cross-sections well away from the fuselage. Moreover the 
lower pressure on the starboard side of the stabilizer suggests that a starboard 
s| de force is imparted to it, thereby creating moments which yaw the aircraft to 
port and roll it to starboard. 

Fig 24 contains a plot of C p distributions in planes EE and GG through the port 
and starboard horizontal stabilizers (cf. Figs. 17 and 19). The effect of the open 
cavity is small on the pressure distribution on the lower surface, but pronounced 
on that of the upper surface, particularly near mid-chord. The C p distributions in 
planes DD and FF of Figs. 17 and 19 are plotted in Fig. 25. The asymmetry of 
the pressure coefficients on the port and starboard stabilizers is less 
pronounced in these cross-sections than on sections EE and GG farther 


inboard. Pressures on the upper surfaces of planes FF and GG on the starboard 
side of the aircraft are lower than pressures on planes DD and EE on the port 
side. This result suggests that a port roll is imparted to the aircraft when the 
telescope cavity is open. 

Cavity and Telescope Sound Pressure Levels 

Fiqs. 26 through 44 contain plots of sound pressure spectra on empennage 
surfaces and in the telescope cavity. The figures depict power spectral density 
(psd) of the sound pressure level (in decibels), defined by 

/ ,2 

spl(db) = 20log,„ 


as a function of frequency, where is the root-mean-square (rms) static 
pressure, in N/m 2 . 

The spectra were computed from 1638 points stored at approximately 6.8 x 10 4 
second intervals in the 1.1 -second data record. The record was zero-padded to 
a lenqth of 2048 points and a Welsh window applied 15 . The Nyquist frequency 
was 739 hz. Hann, Parzen, and Bartlett windowing functions were also applied 
to the data in order to determine what effect, if any, different windows had on the 
spectral representations of the data. The different functions affected some 
spectra at frequencies greater than 100 hz. However, the energy content in 
sound pressure fluctuations at these frequencies is typically several orders of 
magnitude smaller than the energy content of the lower frequency fluctuations 
of primary interest in the design of the SOFIA telescope. Moreover, in t e 
frequency range where the different windowing functions had some effect on 
spectra, numerical dissipation is significant (Atwood 3 ). 

Spectra at points PHI, PH2, PH3, and PH4 on the upper surface of the port 
horizontal stabilizer (of. Figs. 17 and 19) are shown in , Fi 9- 2|- »“ und 
pressure level is greatest at PH4 and least at PH3, but greater at PHI than PH2, 
as is evident from examining the values tabulated in Table 2. This result is 
surprising. One might expect the sound pressure level to be greatest at PH3, 
which is farther upstream and closer to the cavity than the other points. 
However, PH4 is located more directly downstream of the cavity and PH3 is 
shielded from sound emanating from sources deep in the cavity, unlike PHI 
and PH2. 

Spectra at points SHI, SH2, SH3, and SH4 on the upper surface of the 
starboard horizontal stabilizer are shown in Fig. 27. The sound pressure level is 
greatest at SH4 and least at SH3 (cf. Table 2). Not surprisingly, the sound 
pressure levels are greater on the port stabilizer than on the starboard 
stabilizer This result is obvious from examining Table 2 and is made clear by 
replotting the spectra in Figs. 26 and 27 in pairs which represent mirror images 
of one another with respect to the centerline of the fuselage. Spectra at points 
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PHI and SHI are replotted in Fig. 28, at points PH2 and SH2 in Fig. 29, at 
points PH3 and SH3 in Fig. 30, and at points PH4 and SH4 in Fig. 31 . 

Spectra of sound pressure levels on the vertical stabilizer are plotted in Figs 32 
through 35. As is evident from the values listed in Table 2, sound pressure 
levels on the port side of the stabilizer are highest at VP4 and lowest at VP1 
and higher at VP3 than VP2. Point VP3 is the least shielded of all the points 
from sound sources deep in the cavity and the sound pressure level at this point 
is higher than at VP2 and VP3. However, VP4 is closer than VP3 to the 
unsteady wake of the cavity. Not surprisingly, sound pressure levels are higher 
on the port side of the stabilizer than at image points on the starboard side. This 
result is evident from examining the spectra and values tabulated in Table 2. 

Figs. 16 and 20 through 22 contain sketches of the cavity environment and 
locations of points where spectra of sound pressure levels have been 
computed. (Plane HH in Fig. 20 is depicted in Fig. 21.) The spectra are plotted 
in Figs. 36 through 44 and Table 3 contains a tabulation of the sound pressure 
levels at the various points. 

The sound pressure levels at points T2 and T3 on the primary mirror of the 
telescope are identical and the spectra in Fig. 36 are similar. Although 
fluctuating pressures are responsible for structural loads on surfaces, as well as 
sound pressure levels, the similar spectra and identical sound pressure levels 
at these points do not necessarily suggest that there is negligible asymmetric 
loading of the primary mirror. A phase relationship between fluctuations at these 
points would suggest asymmetry in the loading. 

Points SI and T1 are located on the secondary and tertiary mirrors of the 
telescope, respectively. Spectra of sound pressure levels at these points are 
plotted in Fig. 37. The sound pressure level on the tertiary mirror is slightly 
higher than on the secondary mirror, even though T1 is shielded from sound 
sources on the cavity aft ramp in the vicinity of points C4 and A3. Wave 
reflections in the azimuth tube of the instrument may be responsible for this 


Spectra at T1 and T4 and shown in Fig. 38. The sound pressure level on the 
base of the telescope, which is shielded from sound sources in the cavity-port 
shear layer is identical to that at points T2 and T3 on the primary mirror fcf 
Table 2). 7 v ' 

Spectra of sound pressure levels at points C3 and C4 on the aft ramp of the 
cavity are plotted in Fig. 39. The sound pressure level is higher at C3 and the 
power spectrum at that point indicates a significant contribution from a broad 
frequency band with peak at approximately 80 hz. Fig. 40 contains a plot of 
spectra at the center of the aft ramp (point C4). Spectra from the simulation 
reported by Atwood 3 and from the wind-tunnel experiments conducted at NASA 
Ames in 1990 are also plotted. Atwood scaled the experimental data, which 
were taken with the telescope bay mounted forward of the wings, to the flight 
conditions of the CFD simulations. However, neither the wind-tunnel model nor 


the geometry in Atwood's simulation contained the aircraft empennage. A 
frequency peak at 110 hz exists in the experimental data and Atwood's 
simulation; the peak is not evident when the empennage is included in the 
geometry. 

Spectra of sound pressure levels at points A1 , A2, A3, A4, and A5 in Fig. 16 are 
contained in Figs. 41 through 43. Point A1 is located on the fuselage slightly 
upstream of the edge of the telescope cavity. Point A4 is located on the fuselage 
and downstream of the cavity. Points A2, A3, and A5 are located on the cavity 
rim. The sound pressure level at point A1 is, not surprisingly, lower than at the 
other points on the fuselage and rim of the cavity (cf. Table 3). The sound 
pressure levels and spectra at points A3 and A4 (Fig. 42) are similar and exhibit 
the broad-band structure evident in the spectrum at point C3 (Fig. 39) in the 
vicinity of 80 hz. The spectra at points A2 and A5 are plotted in Fig. 43. The 
spectrum at A5 also exhibits the broad-band peak at 80 hz. However, the 
spectrum at A2 contains very energetic peaks at 100 and 200 hz, the latter of 
which is a harmonic of the first. A significant source of sound thus exists on the 
edge of the cavity where A2 is located. As shown in Fig. 20, the edge of the 
cavity is sharp in that region and the sound spectrum at A2 suggests that this 
edge is in an impingement region, unlike the opposite edge in the vicinity of A5, 
which is characterized by a markedly different spectrum. Further analysis 
should focus on determining flow patterns over the fuselage in the vicinity of the 
cavity when it is open to the atmosphere. Doing so may reveal whether an open 
cavity induces upwash over the fuselage behind the wings and the manner in 
which the geometry of the cavity edge in the vicinity of point A2 should be 
modified in order to eliminate the strong source of sound there. 

Examination of the spectra in Fig. 44 reveals that the sound pressure level at 
point C2 on the cavity aft bulkhead is higher than at Cl on the forward 
bulkhead. Moreover, the mean pressure is higher at C2 than at Cl , as indicated 
by the time history of normalized pressure in Fig. 45. This result has an 
important implication for loads on the telescope, as discussed in the next 
subsection. 

Telescope Loads 

Figs. 46 and 47 contain plots of the time histories of integrated force and 
moment coefficients, respectively, on the telescope and Table 4 contains a 
summary of the mean and root-mean-square loads (forces and moments) 
computed from the data record. The lift, drag, and side forces on the telescope 
are components of force in the aircraft (x,y,z) frame of reference. Moments and 
moment coefficients are referenced with respect to the (x'", y'", z") coordinate 
axes, which are parallel to the aircraft axes, but which pass through the center 
of the telescope air bearing. SOFIA design engineers reference telescope 
moments with respect to axes through the air bearing, which is a counter-weight 
mechanism mounted forward of the instrument in the aircraft. The bearing 
isolates the telescope bay from the pressurized sections of the aircraft and must 
support all loads on the telescope, including the weight of the instrument when 
the aircraft is on the ground. The center of the bearing is located on the y = 0 



symmetry plane, 10" forward of the forward bulkhead of the cavity and 6" above 
the longitudinal axis of the fuselage. 

The forces tabulated in Table 4 are small and drag is in the upstream direction 
This result is consistent with the fact that pressure is higher on the aft bulkhead 
of the cavity than on the forward bulkhead (Fig. 45). Power spectra of fluctuating 
forces are contained in Fig. 48. In all spectra a peak containing significant 
energy occurs between 25 and 30 hz, approximately. A peak with much less 
energy occurs near 80 hz. Spectra of moments are plotted in Fig. 49. Significant 
energy occurs in spectral peaks of the pitch and yaw fluctuating moments near 


Open-Cavity Simulation, Wind-Tunnel Conditions 

Results of this simulation were summarized by Machak et al. 16 , a copy of which 
is contained in Appendix B. 


Conclusion 

^report contains a brief history of the computational fluid dynamics effort in 
the SOFIA project and the results of three CFD simulations. Comparisons with 
experiment demonstrate that CFD reproduces well the C p distribution on the 
aircraft fuselage in the vicinity of the empennage, where effects of an open 
telescope cavity on aircraft fatigue and control are likely to be pronounced 
Comparisons with measurements from infrared images of the Shuttle Carrier 
engine exhaust plumes demonstrate that CFD is capable of providing good 
estimates of plume size near the engine exits. CFD provided a better estimate of 
the size of the plume than did a simple model of plume dynamics constructed 
froni empirical correlations. Nevertheless, the computational resources required 
by the latter model are very small and it can be used to generate approximate 
solutions when approximations are deemed sufficient. For example, the model 
might be used to calculate plume temperatures in stray-light analyses. 

The open-cavity simulations provided estimates of sound pressure levels on 
empennage surfaces and in the open telescope cavity. The presence of a 
source of significant acoustic energy on the cavity rim suggests that the cavity 
aperture should be redesigned in the vicinity of the source. The simulations also 
provided estimates of telescope loads and moments. Aerodynamic forces on 
the telescope are small and drag is in the upstream direction. This result is 
consistent with pressure levels in the cavity. At flight conditions significant 
energy exists near 30 hz in all spectra of the fluctuating forces and in the 
spectra of fluctuating yaw and pitch moments. 

The work described in this report has supported a national program of scientific 
importance. As has been discussed, there are a number of technical issues 
which need to be addressed in the SOFIA project and for which CFD has 
proven to be an important design tool. Indeed, since wind-tunnel experiments 
are time-consuming to prepare for and expensive to run, CFD will likely remain 



the pre-eminent design tool available to engineers during the lifetime of the 
SOFIA program. 
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Table 1 Tail Cone Surface Temperatures 


event 

exhaust gas static 
temDerature (°C) 

adiabatic wall 
temperature (°C) 

tail cone temperature 
(°C)' 

i 5 

272 

336 

290 

11 

266 

328 

285 

5 3.6,9,12 

322 

421 

357 
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Table 2 Empennage Sound Pressure Levels 


— 

point 

sound pressure level fdbi 

vemcai staoiirzer ' 1 


VP1 

110. 


VP2 

122. 


VP3 

129. 


VP4 

131. 


VS1 

99. 


VS2 

103. 


VS3 

112. 


VS4 

122. 

port honzontai stabilizer 


PHI 

127. 


PH2 

125. 


PH3 

120. 


PH4 

134. 

starboard horizontal stabilizer 


SHI 

118. 


SH2 

120. 


SH3 

109. 


SH4 

128. 
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Tabie 3 Telescope Cavity Sound Pressure Levels 



point 

sound pressure level (db) 

cavity wails 

Cl 

111 . 


C2 

124. 

cavity aft ramp 

C3 

141. 


C4 

129. 

telescope primary mirror 

T2 

116. 


T3 

116. 

telescope secondary mirror 

SI 

111 . 

telescope tertiary mirror 

T1 

114. 

telescope base 

T4 

116. 

cavity aperture 

A1 

114. 


A2 

140. 


A3 

135. 


A4 

129. 


A5 

133. 
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* F x > Fy, F z drag, side force, and lift in aircraft coordinates (x,y,z) 

Mo^Mo^ Mo roil, pitch, and yaw moments about the center of the 
telescope air bearing; the (x J " ,y"' ,zf") coordinates are 
obtained from the (x,y,z) coordinate system by 
translation 
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Fig. 2 Engine exhaust plume. 







plume width; IR data at 35,000 feet and CFD prediction at 




; empirical and CFD predictions. 






elocity distributions. 








Fig. 9 Plume velocity dis 




Fig. 10 Plume temperature 





Fig. 1 1 Plume velocity distribu 




Fig. 12 Plume temperature distributions at x'/r 




Fig. 13 Plume velocity dislribu 




Fig. 1^1 Plume temperature dis 





Fig. 15 Sketch of the locations of cross-sectional planes AA, BB, and CC 
through the cavity and vertical stabilizer. 
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Fig 1 „ Locations of cross-sections through the horizontal stabilizers and poin 
on their surfaces: (a) starboard stabilizer cross-sections FF and GG and points 

SHI through SH4; (b) port stabilizer cross-sections DD and EE and points PHI 
through PH4. 
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Fig. 20 Cross-section AA through the cavity. 



Fig. 21 Cross-section HH through the cavity. 
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Fig. 23 Vertical stabilizer C p distributions 





Fig. 24 Horizontal stabilizer C distributions on cross-sections EE and GG 
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Fig. 25 Horizontal stabilizer C distributions 




Fig. 26 Spectra of sound pressure levels on the port horizontal stabilizer 
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Fig. 28 Spectra of sound pressure levels at points PHI and $H1 on the 
horizontal stabilizers. 
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Fig. 30 Spectra of sound pressure levels at points PH3 and SH3 on the 
horizontal stabilizers. 







sound pressure level 



Fig. 31 Spectra of sound pressure levels at points PH4 and SH4 on the 
horizontal stabilizers. 





Fig. 32 Spectra of sound pressure levels at points VP1 and VS1 on the vertical 
stabilizer. 
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Fig. 36 Spectra of sound pressure levels at points T2 and T3 on the primary 
mirror of the telescope. 
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Fig. 40 Measured and computed spectra of sound pressure levels at point C4 
on the cavity aft ramp. 
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Fig. 42 Spectra sound pressure levels at points A3 and A4 on. the fuselage 
downstream of the cavity. 
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Fig. 44 Spectra of sound pressure levels at points Cl and C2 on the cavity 
forward and aft bulkheads, respectively. 
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Fig. 48 Power spectra of telescope fluctuating forces. 
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Fig. 49 Power speclra of telescope fluctuating 
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Appendix A 


Empirical Equations Governing the Growth of Coaxial Jet Plumes 

Potential core lengths (cf. Fig. 2): 


*cp / 2r c — 


(4.2- 


UM/tefc 


0.2 


1-0.924 


-2.44 


4W2 

1-0.924 


-2-44 


where M c and M f are the Mach numbers of the hot engine and 
turbofan exhausts, respectively, and 


he =Vf /V c 

ctfe —Tf/T c 
^ if =14 / Vf 
a if I T f 

Plume centerline temperature and velocity ratios: 

T ct 1T C = 1 - exp[l.4 / ( 1-x'/ x' cp )] 
v d I V 0 = 1- exp[l.4 / (1- x' / x' cp )j 

Plume spread rates: 


5 C = Vtc / <j c 
5 f = Vrc / a f 

^ sej = hit / a se j 

5f/j = 1-46 se y 

where S th is the spread rate of the thermal plume in a single-equivalent jet 


A-1 


Pc = tan 1 5 C 
P f = tan -1 5 f 
Pse/ = tan 5 sei 

Ptf» = tan -1 5 th 


where, if X fc < 0.8, 


10.0 


<*c = 


(1-0.29/W c 2 )[l + 0.5(/W c 2 -1)(7 C !T f -14) 2 ] (1-0.924) 

10.0 

Gf = (1-0.29/W f 2 )[l + 0.5(M f 2 -X)(T f /T„ -14) 2 ] (1-0.92^) 


10.0 


®sej — 


(1-0.29/W c 2 )[l + 0.5(M c 2 -1)(7 c /T„ -1.4) 2 ](1- 0.924) 


or, if X fc > 0.8, 


10.0 


°c = 


<j f = 


&sej ~ 


(1-0.29/W c 2 )[l + 0.5(M c 2 -1)(T C /Tf -1.4) 2 ] 

mo 

(1-0.29A4 2 )[l + 0.5(/W f 2 -1)(7, /T. -U) 2 J 

mo 

(1 - 0.29 M 2 )[l + 0.5(M C 2 - 1)(7 C I T - - W 


Single-equivalent jet radius: 


r sej — - n + 


r f 2 T c V f 2 


r c 2 T f V 2 


Radial velocity and temperature distributions (fully-developed jet): 

V ~V„ 


Vc ,~ V - 


7-7, 
T c i — 7„ 


= 0.9 


1- sin-r-j 


7T77 I 

TJ' 


I ;TT7 1 

= 0.9J- sin— j , 


- 1 < 77 <1 


-1<77 <1 
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where 


Root-mean 


where 

h 


T l=(z'-z '6.5) /zas 

square temperature distribution (fully-developed jet): 
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Abstract 

A survey is done of previous and current testing and 
analysis on the aerodynamic configuration of the 
Stratospheric Observatory for Infrared Astronomy 
(SOFIA). Current studies use both wind tunnel testing and 
computational fluid dynamics (CFD) to determine the 
SOFIA baseline cavity/shear layer control design. A 
comparison of candidate shear layer control schemes is 
performed for the open cavity. These schemes are 
evaluated on the basis of cavity acoustic levels, impact to 
flow over the tail control surfaces, amount of 
aerodynamically induced torques on the telescope, and 
aero-optical quality of the shear layer flow. A comparison 
of the 747-SP and 747-200 as a suitable platform for 
SOFIA is also addressed using the above criteria. 

Background 

SOFIA is a planned 2.5 meter aperture Cassegrain 
telescope that will be installed in a Boeing 747 aircraft 
and operated at altitudes from 41,000 to 46,000 feet. It is a 
follow-on mission to NASA's Kuiper Airborne 
Observatory (KAO), which consists of a 0.9 meter 
telescope in an open-port (no material window) cavity in a 
Lockheed C-141. An airborne observatory offers a 
number of advantages over both ground-based and space- 
based telescopes. Much of the infrared spectrum is 
absorbed by water vapor in the atmosphere, but SOFIA 
will operate above nearly all atmospheric water vapor, so 
a much broader spectrum of wavelengths is available for 
observation. Another advantage of an airborne platform is 
the ability to take advantage of celestial events of 
opportunity, such as novae and planetary occultations, 
which may not be accessible by ground based 
observatories. Also, compared to space-based 
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observatories, the SOFIA mission can be carried out at a 
fraction of the cost. However, an airborne observatory 
creates some challenging technical issues that need to be 
addressed. 


Figure 1: The SOFIA Aircraft 

To diminish IR transmission losses, SOFIA will operate 
with an open-port cavity. This creates a set of challenging 
aerodynamic/aeroacoustic design problems. The shear 
layer over the cavity aperture must be controlled to 
eliminate unwanted cavity resonance and reduce the 
overall cavity sound pressure level. This must be done in 
a way which impacts the telescope's optical and 
mechanical performance as little as possible. 

Effects of an open cavity 

There are several specific areas where the performance of 
the SOFIA aircraft is affected by the presence of an 
open-port cavity. These areas are discussed here in 
general terms before addressing how these challenges are 
being met through a combination of testing and analysis. 

Cavity Resonance 

Rectangular cavities will tend to resonate when excited by 
air flow. From a safety and reliability standpoint, an 
aircraft operating with an open cavity must incorporate 
some method of eliminating this resonance. A resonant 
cavity could affect the stability of the aircraft and increase 
the fatigue of components in and near the cavity. An 
acoustic resonance would also increase the level of force 
on the telescope. A primary goal of the wind tunnel tests 
and CFD analysis is to identify shear layer control 
schemes which provide the lowest possible sound 
pressure level (SPL) in the cavity. 
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Telescope loads 

Unsteady pressure fluctuations in the cavity have the 
capacity to produce forces on the telescope which must 
then be counter-acted by the torque motors. As the 
amplitude of the aerodynamically induced torque is 
increased, the size (and therefore the cost) of the torque 
motors must increase accordingly. These aerodynamic 
forces will also impact the pointing accuracy of the 
telescope. _ 

Aircraft Stability and Control 

Originally, the cavity was located in the forward part of 
the airplane, just in front of the wings. With this concept, 
the stability and control of the aircraft was not a concern, 
since the shear layer would have almost the entire length 
of the aircraft to re -attach and stabilize before 
encountering the control surfaces on the empennage. The 
cavity location was subsequently moved to the back of the 
aircraft as a cost saving measure. Now that the cavity is 
located aft of the wings, a greater potential exists that the 
flow over the tail surfaces will be altered from that of a 
clean configuration (no-cavity). 

Aircraft Drag 

The cavity and shear layer control device add drag to the 
aircraft. This is a concern because the drag makes a direct 
impact on the available time at 41,000, which is a top 
level requirement. Estimates from contractor studies 
predict that the coefficient of drag (Co) will increase by 
.0015 due to the cavity and shear layer. This would result 
in the time at flight being reduced by approximately 12%, 
or 1.5 hours. Currently, there appears to be a considerable 
margin for the time at altitude requirement, but these 
estimates have not been verified by either experiment or 
by computation. Even if the minimum time at altitude 
requirement is met, changes in drag could have an effect 
on the operational costs of the observatory. 

Seeing 

An open port cavity impacts the quality of the image. As 
the high speed flow passes over the cavity, the transition 
region from freestream condition to the cavity conditions 
is referred to as the shear layer. The shear layer is highly 
turbulent which means that there will be a high level of 
density fluctuations present. The index of refraction of air 
is related to the density, so as the light passes through the 
shear layer, it is getting refracted by the density 
fluctuations. Temperature gradients and 'hot spots' within 
the cavity can create an additional error. Both the shear 


layer and cavity conditions result in a wave front error of 
the image when it reaches the focal plane. 

1990 wind tunnel test 

A wind tunnel test was performed in 1990 with the main 
objective being to determine whether flow control devices 
can suppress cavity resonance. Secondary objectives were 
to optimize the shear layer control device to give the 
lowest cavity sound pressure level and to determine the 
aero-optical quality of the flow in the shear layer. All 
these objectives were met during the test. The model used 
for this test was a 7% scale model of a 747-SP (Fig. 2) 
which was sting mounted and had the wings clipped just 
beyond t he inboard eng ines. 



Figure 2: Forward Cavity Wind Tunnel Model 
747-SP, 1990 Test 


Test Conditions 

• Mach number (Moo): 0.63 to 0,88 (0.85 reference) 

• Reynolds number/foot: 3.8 x 10^ 

• Stagnation conditions: 14.7 psia, 580° R 

• Angle of attack(a): 0° to 3° (2.5° reference) 

• Yaw angle(P): ±1° (0° reference) 

Results and Discussion 

The 1990 test resulted in a shear layer control concept 
which had a quasi-elliptical aperture shape, with a molded 
ramp on the downstream edge of the cavity * . This 
solution was arrived at after considering other treatments 
such as straight aft ramps, various devices such as 
upstream porous fences, and combinations of both. The 
final quieted cavity configuration had a sound pressure 
level(SPL) of 137.2 dB wind tunnel scale, which would 
translate to approximately 127 dB full scale. Aero-optical 
measurements showed that the molded quasi-elliptical 
configuration performed much better for seeing than using 
a porous fence upstream of the cavity, particularly for 
wavelengths less than 10 \im. Above 10(im, the telescope 
is diffraction limited and shear layer seeing is not a factor. 
The porous fence is used by the Kuiper Airborne 
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Figure 3: 747-200 Wind Tunnel Test Model 
from 1 994 Test 

Observatory (KAO) to thicken the boundary layer and 
eliminate the cavity resonance. 

1994 wind tunnel test 

During the summer of 1994, a second test was held. The 
Principal Investigator for this test, as in the 1990 test, was 
Dr. William Rose, of Rose Engineering and Research, Inc. 
The test was carried out using the same test conditions as 
the 1990 test and a modified test modeK 2 ' 3 ). The model 
from the previous test was modified to add the 
empennage, move the cavity location to the aft of the 
aircraft, and to extend the wings to the location of the 
outboard engines. The new model can be configured as 
either a 747-SP or a 747-200. Fig. 3 shows the model 
configured as a 747-200, and Fig. 4 shows the model 
configured as a 747-SP. A remotely controlled telescope 
and cavity door provided for multiple telescope positions 
during each run. The test conditions were the same as 
during the 1990 test. 

Instrumentation 

The following instrumentation was used during the test. 

• Pitot and static pressure ports for steady pressure 
measurements 

• Kulite pressure transducers for unsteady pressure 
measurements 

• Moveable rake drive for pitot probes in the shear layer 

• Flow visualization: tufts, oil paint flow 

• A 2.5" 6-component balance to measure tail loads 

• A 0.75" 6-component balance to measure telescope 
forces and moments 

• Accelerometers mounted on the telescope used for 
measurement of aerodynamic loads on the telescope. 

• Unsteady pressure transducers mounted on the 
telescope to measure higher frequency pressure 


Figure 4: 747-SP Wind Tunnel Test Model 
from 1994 Test 


variations which can be analyzed to predict forces and 
moments. 

• Pressure Sensitive Paint (PSP) applied to the port side 
of the vertical stabilizer and to the top of the port 
horizontal stabilizer, to examine the effect of the open 
cavity on tail pressures. 

Telescope Load Measurement 

The telescope loads were measured by two independent 
methods, both using instrumentation mounted on the 
telescope model, shown in Fig. 5. The 'pressure method' 
derived the telescope loads from a measurement of the 
unsteady pressures on the surface of the telescope. The 
'dynamic model method’ 4 calculated the force being 
applied to the telescope by modeling the dynamic 
response of the telescope and using outputs from 



Figure 5: Wind Tunnel Telescope Model 


3 

American Institute of Aeronautics and Astronautics 





origin +Z - Telescope Line»Of 

Sight (LOS) Axis 4 


+X - Telescope Elevation Axis, 

+Y - (out of page) Telescope 
Cross-Elevation Axis 

Figure 6: Telescope Assembly Coordinate System 
(Fixed to Telescope Frame of Reference) 

accelerometers mounted on the telescope and a 
6-component balance upon which the telescope was 
mounted. Both methods only measured the loads acting in 
the direction of the telescope Line of Sight (LOS) axis, 
which produce moments about the Cross-Elevation(CE) 
axis. This coordinate system is shown in Fig. 6. Due to the 
fact that this rotation axis has the largest surface area 
capable of generating a moment, it is expected that 
aerodynamically induced moments about this axis will the 
highest. 

The pressure measurement method used the output from 8 
unsteady pressure transducers mounted on the telescope, 4 
on the top of the primary mirror surface, and 4 on the 
bottom. The data was collected digitally and downloaded 
to a PC for processing. From there, an average pressure 
difference was calculated at each time step to determine a 
torque about the CE axis. 

The dynamic measurement method was implemented by 
Paul Keas of SVERDRUP. Two accelerometers were 
mounted on the telescope which could measure 
accelerations along the LOS axis. Also measured were the 
forces imparted to the telescope from the bulkhead by 
means of a 6-component balance. Newton's second law, 
the sum of the forces equals mass times acceleration, was 


used to determine the force resulting from fluctuating 
pressures. An assumption is made that the only forces 
acting on the telescope are from the acoustics and the 
reaction forces from the balance. The telescope mass is 
known, and the balance force and telescope acceleration 
are measured. From there it is possible to calculate the 
acoustic force. 


Results and Discussion 
Apemire^comnarison 

Four basic shear layer control concepts were tested on the 
747-200, plus one case where there was no treatment of 
the aperture. These treatments are shown in Figures 8-12. 
In addition to the rectangular and D-shape apertures, a 
splitter plate aperture was tested on the 747-SP (Fig, 13). 
The configuration numbers refer to the order in which the 
different design concepts were tested. Each configuration 
has a unique aperture and cavity treatment. 

Although the circular aperture (Fig. 9) design had the best 
performance on the forward cavity configuration, the 
same was not true in the aft of the aircraft. The upstream 
curved section had more of a benefit in the forward cavity 
case due to the presence of non-parallel, angular flow 
caused by the proximity of the wings. In the aft location, 
the flow is mostly parallel and has only a slight 
angularity. An improvement in the cavity SPL was 
obtained by making the upstream portion of the aperture 
straight while keeping the downstream section curved. 
This resulted in the D-shape aperture, shown in Fig. 12. 
Configuration 65, shown here, was the best solution tested 
for the 747-200. The same aperture shape was 
subsequently tested on the 747-SP platform in 
configuration 86 and produced similar results. Having a 
curved downstream edge means that the aperture will 
have to move with the telescope through the 20° - 60° 
elevation range. To examine less costly options, two 
straight ramp designs were investigated. In both of these 
designs, the ramp would remain fixed to the aircraft, and 
the door design would be simpler. However, both the 
internal straight ramp (Fig. 11) and the external straight 
ramp (Fig. 10) had unacceptable performance. The splitter 
plate configuration was designed to create a thinner shear 
layer over the cavity and therefore improve seeing. 
However, due to higher cavity sound pressure level and an 
adverse effect on flow over the tail, this design was 
deemed unacceptable. 

Fig. 14 shows the overall cavity SPL for each of the 
different types of aperture shapes tested. 3 Fig. 15 shows 
the power spectral density function of the pressure on the 
forward bulkhead for configuration 86. As we can see 
from Fig. 14, a 10 dB reduction in the cavity SPL is 
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obtained when comparing the D-shape aperture to the 
untreated cavity. For the better configurations (i.e. 65, 86) 
the cavity SPL varied little with respect to either Mach 
number, angle of attack, yaw angle, or telescope elevation 
angle. This shows that the design is robust and should 
maintain its performance through changing conditions. 

Telescope Torque M easurement 

The loads generated on the telescope were calculated for 
all configurations tested during the wind tunnel test using 
the pressure method. The dynamic method could not be 
used for some of the aperture treatments due to a higher 
risk of damage to the balance. The telescope balance was 
locked out on these configurations. Results from the 
pressure method are presented in Figures 16-19. 

The induced moments about the telescope cross-elevation 
axis are plotted against the aircraft angle of attack in 
Fig. 16 for a variety of different aperture treatments. In 
addition to the six configurations shown in Figures 8-13, 
tfiere is also configuration 86, which is the D-shape cavity 
implemented on the 747-SP. As can be seen from the 
graph, the D-shape aperture configurations not only have 
the lowest RMS moments, but also exhibit virtually no 
dependence on angle of attack. The plot for the untreated 
aperture (configuration 18) demonstrates the tremendous 
improvements which can be obtained using passive shear 
layer treatments. 

Fig. 17 shows the cross-elevation moments plotted as a 
function of Mach number for the same configurations. 
The general trend is that the moments will tend to 
decrease at higher Mach numbers, especially for the 
D-shape apertures. Fig. 18 shows the dependence on yaw 
angle of the aircraft for the D-shape apertures on the 
747-200 and 747-SP. At cruise conditions, there is 
virtually no effect of a yaw angle on telescope cross- 
elevation moments. 

Figures 16-18 show that there is little variation in the 
telescope cross-elevation moments throughout a wide 
range of conditions. For the better configurations, there 
was also very little influence on the moments due to 
changes in telescope elevation angle. 

The cross-elevation moments are plotted against the 
overall cavity sound pressure level in Fig. 19. Each point 
represents one configuration at cruise conditions. Only 
configurations which were variations on the D-shape 
aperture are plotted here. The dashed line in the figure 
represents the best fit to the data assuming a linear trend. 
The slope of this line is approximately 91 ft-lb/dB, 
However, as can be seen at the low end of the curve, the 
slope may be starting to flatten out. However, it appears 


that decreasing the cavity SPL even 1-2 dB may have a 
significant impact on the telescope moments. 

Figures 20 and 21 present a comparison of the two 
methods used for measuring the cross-elevation moments. 
Fig. 20 shows a PSD of the cross-elevation moments and 
Fig. 21 also shows a PSD for the same wind tunnel test 
sequence. Note that both methods predict a large spike in 
the spectrum right at 400 Hz. Smaller spikes occur at 
around 250 Hz, and around 900 Hz. The spike in the 
spectrum at 35 Hz in the dynamic model spectrum is 
attributable to the telescope balance natural frequency. 
Given that both methods of measurement are completely 
independent of one another, this leads to a high degree of 
confidence in the results. 

To get an idea of why a spike occurs at 400 Hz in the CE 
moment spectrum, compare Figures 15 and 24, which 
show pressure PSDs on the forward bulkhead and the 
primary mirror respectively. Both plots are for the same 
configuration and flight conditions. The forward bulkhead 
spectrum has energy at many frequencies. However, the 
structure of the telescope acts in a way to preserve the 
energy around 400 Hz, and reduce the acoustic energy at 
other frequencies. This leads to the postulation that an 
open, truss structure telescope may act to reduce the 
acoustic loads. This will be investigated during a wind 
tunnel test scheduled for January 1995. 

Aircraft St ability and Control 

The issue of aircraft stability and control was addressed 
during the 1994 wind tunnel test. This was done by 
measuring the loads on the empennage and the pressure 
distribution over the tail surfaces. The Principal 
Investigator for this portion of the test was Dr. Ian 
Gilchrist of Analytical Methods Inc. Preliminary analysis 
of the results indicate that all the internal aperture 
treatments had a negligible effect on tail lift. 5 Effects of 
the cavity on directional and longitudinal stability will be 
negligible at cruise conditions. 5 Stability of the aircraft at 
lower speeds has not yet been addressed by the project. A 
plan on how to deal with this issue is being developed. 

Aircraft C omparison 

From the wind tunnel tests, it appears that the presence of 
an open-port cavity had the same effect on both the 
747-200 and 747-SP platforms. This is true when 
comparing the effect that similar aperture treatments had 
on cavity acoustics 3 , telescope moments, and aircraft 
stability and control 5 . The decision on which aircraft to 
use may then be based on questions of design complexity 
and operational costs for the two aircraft. 
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Computational Fluid Dynamics (CFD) 
Background 

The investigation of the aerodynamic issues pertinent to 
the SOFIA design has progressed along parallel paths 
using wind tunnel testing and computational simula tio ns. 
The science of computational fluid dynamics (CFD) has 
developed to the point where very complex flows, both 
external and internal, can be simulated to an accuracy 
comparable to the experimental measurements. It is 
desirable to have good CFD techniques to complement 
experiments, especially in areas where precise 
experimental measurements are difficult. The SOFIA 
cavity is one such environment, where the flowfield is 
complicated by a complex aircraft platform and a 
resonating cavity housing and a telescope located below 
an unsteady turbulent shear layer. Therefore, the twofold 
objective of the numerical simulation of the flow about 
SOFIA airborne observatory are 1) to complement the 
wind tunnel experiments in the safe design of not only the 
aircraft platform that carries the telescope but also to have 
the least detrimental effects upon the telescope optical and 
mechanical performance; and 2) to provide the detailed 
flowfield information that are difficult to measure in the 
wind tunnel. 

The SOFIA CFD program started in 1989. At first, simple 
2-D geometries were considered for the cavity and 
modeled to verify the numerical schemes that would be 
used in the full-scale problem. During the first year, work 
began on simulation of the flowfield over the 747-SP 
geometry, the cavity, and the telescope for the forward 
cavity configuration. By the end of 1991, flowfield 
solutions were obtained for a sting-mounted 747-SP with 
and without a forward-located (with respect to the wing) 
cavity 6 * 7 . This was the same geometry that was used in 
the 1990 wind tunnel experiment. Calculations were 
performed for two different cavity cases and the results 
compared favorably with wind tunnel experiments. In 
1992, the model geometry was changed to reflect the re- 
location of the cavity to the aft (downstream of wing) end 
of the aircraft®' 10 . These first simulations of the aft cavity 
used a sting mounted (i.e., no tail) configuration, and the 
later calculations included the tail empennage. The current 
CFD work on SOFIA investigates the aft cavity 
configuration for both the 747-SP and 747-200 aircraft. In 
the current simulations, the same circular aperture 
treatment is used that resulted from the 1990 wind tunnel 
test. Future simulations will incorporate the D-shape 
aperture treatment from the 1994 wind tunnel test. 


Method of Solution 

The numerical method solves the Reynolds averaged 
Navier-Stokes equations in the thin layer limit on an 
overset grid framework using the flow solver 
OVERFLOW 11 This flow solver uses a central- 
differenced, implicit, diagonal scheme 12 with added 
second- and fourth-order numerical dissipation terms. The 
numerical scheme is second-order accurate in space and 
first-order accurate in time. The flow on the entire 
configuration is assumed to be fully turbulent at the 
Reynolds number of interest and a simple algebraic 
turbulence model 13 is used to estimate the eddy viscosity. 
A no-slip boundary condition is specified at the wall with 
zero normal pressure gradient along with an adiabatic 
wall condition. To update the information exchange at the 
overset grid interfaces, a trilinear interpolation of the 
dependent variables is used. The computations were 
performed on the NAS Cray C-90 and the ACSF Cray 
C-90 supercomputers at NASA Ames Research Center, 



Figure 7: Representation of Numerical Grids 
on the 747-200 Computational Model 


The gridding of the SOFIA platforms of the Boeing 
747-SP and 747-200 aircraft is a tedious job. Initially, the 
surface grids were reconstructed on several individual 
surfaces from CAD data using the GRIDGEN 14 and 
S3D 15 codes. Once the complete surface (with overlaps) 
grids are constructed to the accuracy needed, then 
structured volume grids are carefully gen erated o ver the 
individual surfaces. Grids of O-O, C-O, and H-H topology 
have been constructed using the HYPGEN grid 
generator 16 . In addition, special tools developed at NASA 
Ames for constructing collar grids and cap grids were also 
utilized. Typically the clean configuration has 38 different 
overset grids and contains over 2.5 million grid points. 
Only engine pylons were not included in this 
configuration. For example. Fig. 7 shows the surface grids 
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for the 747-200 configuration without cavity. The cavity 
comprises of 7 additional grids and contains about 1 .6 
million grid points. After the surface and volume grids 
are constructed, then hole boundaries are created by 
cutti "g holes in adjacent grids using the PEGSUS 
code • . This is a tedious iterative process and requires 
several iterations of hole punching and grid generation to 
eliminate all orphan points in the flowfield. This hole 
boundary information is necessary to run the flow solver. 
For the 747-200 there are a total of 45 grids including the 
cavity grids and they consist of over 4 million grid points. 
After punching holes in the overset grids, the calculations 
were run in a unsteady fashion by starting from the steady 
solution. 

Simulation Flow Conditions 

• Mach number: 0.85 

• Reynolds number: 4.0x1 0 6 /foot 

• Angle of attack: 2.5° 

• Yaw angle: 0° 


Results and Discussion 

Figures 22 - 26 show comparisons of CFD results with 
measurements made in the wind tunnel. For these cases, 
configuration 39 from the wind tunnel is used for 
comparison to the 747-200 solution, and configuration 86 
is used for comparison to the 747-SP simulation. Both 
simulations were performed at wind tunnel conditions on 
models that had the wings clipped at the location of the 
outboard engine, just as the wind tunnel model was 
configured. Note that neither of the wind tunnel 
configurations match exactly the CFD configurations. 
Configuration 39 has a circular type aperture whose 
overall dimensions were closest to those of the CFD 
simulation. Configuration 86 is the final D-shape aperture. 
No circular apertures were tested on the 747-SP in the 
wind tunnel. 

Figures 22 and 23 show power spectral density plots of 
points on the forward bulkhead. The experimental records 


from which the spectra were computed were five seconds 
long. The records from which the CFD spectra were 
computed were approximately .07 seconds long. The level 
of the spectra compares well with the experimental data in 
the 400 Hz to 1000 Hz range for wind tunnel conditions. 
Below that range the CFD simulation is affected by the 
brevity of the simulation and above that range numerical 
dissipation causes a fall-off of the spectrum. In Figures 24 
and 25, the power spectrum is compared for a point on the 
surface of the primary mirror. The overall levels of the 
CFD prediction are below that of the experiment, but the 
major peaks in the spectrum are captured, particularly at 
400 Hz. Fig. 26 contains a comparison table of wind 
tunnel and CFD results for the sound pressure level at 
points on the telescope and bulkhead surfaces. 

Faff 27 presents a plot of instantaneous Mach number 
contours for the 747-SP simulation in a 2-D cross-section 
of the cavity, perpendicular to the cross-elevation axis of 
the telescope. The flow within most of the cavity is below 
Mach 0.1 except for a few locations near the comers of 
the cavity and telescope. Fig. 28 shows a mapping of 
sound pressure levels on the forward and aft bulkheads. 

Conclusion 

A comparison of results from the wind tunnel experiment 
and CFD simulations provides verification of CFD as a 
prediction tool. Wind tunnel tests over various 
configurations and flight conditions have proven the 
feasibility of locating the cavity in the aft of the aircraft 
with considerable margin. The data show that there will 
be no impact to the aircraft stability and control at cruise 
conditions. The wind tunnel test has proven that a passive 
quieted cavity configuration can be built and a baseline 
aperture treatment has been identified. The 
aerodynamically induced moments about the cross- 
elevation axis are larger than desired, but some possible 
solutions to this problem are being investigated. The wind 
tunnel tests and CFD simulations have shown that either 
the 747-SP or the 747-200 are viable candidates for the 
SOFIA platform. 
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Figure 8 Circular Aperture 
Configuration 8 from 1994 Wind Tunnel Test 



Figure 9: Untreated Aperture 
Configuration 18 from 1994 Wind Tunnel Test 



Figure 10: External Ramp 
Configuration 19 from 1994 Wind Tunnel Test 



Figure 1 1 : Straight Internal Ramp 
Configuration 31 from 1994 Wind Tunnel Test 



Figure 12: D-Shape Aperture 
Configuration 65 from 1994 Wind Tunnel Test 



Figure 13: Splitter Plate 
Configuration 92 from 1994 Wind Tunnel Test 
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Figure 14: Comparison of Overall Cavity Sound Pressure Level for Different Configurations 
Moo=.85, ot=2.5°, p=0°, Wind Tunnel Scale 
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Figure 15: Power Spectral Density (PSD) of Forward Bulkhead Pressure for 747 -SP 
Configuration 86, Moo=.85, a=2.5°, 0=0°. Wind Tunnel Scale 
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Figure 16: Effect of Angle of Attack on Telescope Cross-Elevation Moments 
Plotted by Configuration Number, Moo—. 85, P=0°> Scaled to Full Scale Flight Conditions 
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Figure 17: Effect of Mach Number on Telescope Cross-Elevation Moments 
Plotted by Configuration Number, a=2.5°, P=0°t Scaled to Full Scale Flight Conditions 
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Figure 18: Effect of Mach Number on Telescope Cross-Elevation Moments 
Plotted by Configuration Number and Yaw Angle, a=2.5°, Scaled to Full Scale Flight Conditions 
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Figure 19: Effect of Cavity Sound Pressure Level on Telescope Cross-Elevation Moments 
D-shape aperture treatments, Moc=.85, a=2.5°, 3=0°, Scaled to Full Scale Flight Conditions 
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Figure 20: PSD of Cross-Elevation Moments for 747 -SP using 
essure Measurement Method, Configuration 86, Moo=.85, oc=2.5°, p=0°. Data in wind tunnel scale 
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Figure 2 1 : PSD of Cross-Elevation Moments for 747 -SP using 
Dynamic Model Method, Configuration 86, Moo=.85, 0=2.5°, p=0°. Data in wind tunnel scale 
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Figure 22: Aft Bulkhead Power Spectral Density Comparison for the 747 -SP Platform, M«,=.85, a=2.5° 
CFD simulation at wind tunnel conditions, Wind Tunnel Configuration 86 (D-shape Aperture) 
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Figure 23: Aft Bulkhead Power Spectral Density Comparison for the 747-200 Platform, Moo=.85, a=2.5° 
CFD simulation at wind tunnel conditions, Wind Tunnel Configuration 39 (Circular Aperture) 
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24: Telescope Primary Mirror Power Spectral Density Comparison for the 747 -SP Platform, Moo- 85, a-2.5° 
CFD simulation at wind tunnel conditions, Wind Tunnel Configuration 86 (D-shape Aperture) 



Figure 25: Telescope Primary Mirror Power Spectral Density Comparison for the 747-200 Platform, Moo— 85, 0-2.5 
CFD simulation at wind tunnel conditions, Wind Tunnel Configuration 39 (Circular Aperture) 
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Figure 26: Comparison of Sound Pressure Level between CFD and Wind Tunnel 
Configuration 39 used for 747-200 comparison, Configuration 86 used for 747-SP comparison, M«,=.85, a=2.5°, (3=0° 


15 

American Institute of Aeronautics and Astronautics 
































































Figure 27: Instantaneous Mach Number Contours from 747 -SP Simulation at 
Wind Tunnel Conditions, Moo = *85, ot=2.5 , P^)° 
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Figure 28: Sound Pressure Levels on Bulkhead Surfaces, Wind Tunnel Simulations, Moo=.85, a=2.5°, (J=0° 
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